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The single-particle spectral function of 56 Ni has been computed within the framework of self-consistent 
Green's functions theory. The Faddeev random phase approximation method and the G-matrix technique are 
used to account for the effects of long- and short-range physics on the spectral distribution. Large scale calcu- 
lations have been performed in spaces including up to ten oscillator shells. The chiral N 3 LO interaction is used 
together with a monopole correction that accounts for eventual missing three-nucleon forces. The single-particle 
energies associated with nucleon transfer to valence 1 pOf orbits are found to be almost converged with respect 
to both the size of the model space and the oscillator frequency. The results support that 56 Ni is a good doubly 
magic nucleus. The absolute spectroscopic factors to the valence states on A = 55, 57 are also obtained. For 
the transition between the ground states of 57 Ni and 56 Ni, the calculations nicely agree with heavy-ion knockout 
experiments. 

PACS numbers: 31.10.+z,31.15.Ar 



I. INTRODUCTION 



The way shell closures and single-particle energies evolve 
as functions of the number of nucleons is presently one of the 
greatest challenges to our understanding of the basic features 
of nuclei. Doubly-magic nuclei are particularly important and 
closed shell nuclei like 56 Ni and 100 Sn have been the focus of 
several experiments during the last years fU, B HI 0] ■ 
Their structure provides important information on theoret- 
ical interpretations and our basic understandingof matter. 
In particular, recent experiments 01 01U, 0] have 
aimed at extracting information about single-particle degrees 
of freedom in the vicinity of 56 Ni. Experimental information 
from single-nucleon transfer reactions and magnetic moments 
JH 0HL 0] > can be used to extract and interpret compli- 
cated many-body wave functions in terms effective single- 
particle degrees of freedom. Transfer reactions provide for ex- 
ample information about the angular distributions, the excita- 
tion energies and the spectroscopic factors of possible single- 
particles states. If one can infer from experimental data that 
a single-particle picture is a viable starting point for inter- 
preting a closed-shell nucleus like 56 Ni, one can use this nu- 
cleus as a basis for constructing valence-space effective inter- 
actions. These interactions can in turn be used in shell-model 
calculations of nuclei with several valence nucleons above the 
N = 28 and Z = 28 filled shells of 56 Ni. Recent measure- 
ments of spectroscopic factors of 57 Ni in high-energy knock- 
out reactions |0] seem to indicate that low-lying states in 57 Ni 
can be characterized as single-particle states on top of 56 Ni 
as a closed-core nucleus. Large-scale shell-model calcula- 
tions by Horoi et al |0] corroborate these findings, whereas 
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a recent experiment on magnetic moments of the ground state 
of 57 Cu Ja], expected to be described as one valence proton 
outside 56 Ni, resulted in much smaller moments than those 
expected from a single-particle picture. Similarly, large tran- 
sition matrix elements between the 0[ ground state and the 
first excited 2[ state indicate that the 56 Ni core is rather soft 
|0], or stated differently, it implies a rather fragmented single- 
particle picture. On the other hand, one ought keep in mind 
that quenchings of spectroscopic factors to about 60% are 
common even for good closed shell nuclei 10]. Experimen- 
tally, spectroscopic factors are defined as the ratio of the ob- 
served reaction rate with respect to the same rate calculated as- 
suming a full occupation of the relevant single-particle states. 
They are therefore often interpreted as a measure of the oc- 
cupancy of a specific single-particle state. However, from a 
strict theoretical point of view spectroscopic factors are not 
occupation numbers but a measure of what fraction of the full 
wave function can be factorized into a correlated state (often 
chosen to be a given closed-shell core) and an independent 
single-particle or single-hole state. Large deviations from the 
values predicted by an independent-particle model, point to a 
strongly correlated system. In this regime, collective excita- 
tions that behave like single-particle degrees of freedom — that 
is quasiparticles — can still arise. 

The above mentioned large-scale shell-model calculations 
|0] have been performed in one major shell, the 1 p0/-shell, 
with an effective interaction fitted to reproduce properties of 
several nuclei that can be interpreted in terms of these single- 
particle states. The number of possible Slater determinants 
that can be constructed when distributing eight valence pro- 
tons and eight valence nucleons in the IpOf shell is more 
than 10 9 . This means that the inclusion of more compli- 
cated particle-hole excitations from shells below and above 
the 1 pOf shell, are well beyondpresent capabilities of large- 
scale diagonalization methods 
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eludes as many as possible of these neglected particle-hole 
excitations. However, there are other many-body methods 
that allow for a computational scheme which accounts for 
a systematic inclusion of more complicated many-body cor- 
rections. Typical examples of such many-body methods are 
coupled-cluster metho ds | Til fl5i jl6l [l7ll . various types of 
Monte Carlo methods 11181 [l9l 12011 . perturbative many-body 
expansio ns 12 ll 12211. self-consistent Green's functions (SCGF) 
methods jllHSlllSlgll^J^Jll, the density-matrix 
renormalization group (3 
functional theory I Bo, [3 
able methods. 
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and the no-core shell-model approach 111 31 14 11 
l46ll have been successfully applied to the theoretical descrip- 
tion of light nuclei with mass numbers A < 12, and Hamil- 
tonians based on nucleon-nucleon and three-nucleon interac- 
tions. However, present experimental studies of nuclear stabil- 
ities are now being pushed to larger mass regions, with mass 
numbers from A — 40 to A = 100. Traditionally, this has 
been the realm of the nuclear shell-model and nuclear density- 
functional theory. These methods employ Hamiltonians and 
density functionals with phenomenological corrections and 
are not directly related to the vacuum nucleon-nucleon in- 
teraction employed in ab initio calculations (exceptions are 
found when perturbative many-body methods are used i22tl ). 
However, in selected medium-mass nuclei, ab initio structure 
calculations can be performed using approaches like coupled- 
cluster and Green's functions theories. These methods allow 
studying ground- and excited-state properties of systems with 
dimensionalities beyond the capability of present large-scale 
diagonalization approaches, with a much smaller numerical 
effort when compared to diagonalization methods aiming at 
similar accuracies. The accuracy of these methods is suffi- 
ciently high to attribute an eventual disagreement between ex- 
perimental data and theoretical results to missing physics in 
the Hamiltonian. In this way, such calculations help to in- 
crease our understanding of the nuclear interaction on a very 
fundamental level. 

Recent coupled-cluster calculations \17{] have reported 
practically converged results of the ground state of medium- 
mass nuclei like 40 Ca, 48 Ni and 48 Ca using the bare chiral 
interaction N 3 LO ll47ll . These calculations employed a har- 
monic oscillator basis to construct the single-particle basis 
and included correlations of the so-called singles and dou- 
bles types. It means that one-particle-one-hole and two- 
particle-two-hole correlations acting on a many-body Slater 
determinant were summed to infinite order. Recently, results 
with three-particle-three-hole correlations have also been ob- 
tained 14811 . The calculations were performed in a harmonic 
oscillator basis containing up to fifteen major shells and re- 
sulted in basically converged ground state properties for a 
given Hamiltonian. 

Another method with a strong potential for performing ab 
initio calculations of nuclei beyond A = 12 is self-consistent 
Green's functions theory. Differently from the coupled-cluster 
and no-core shell model, this method does not construct the 
wave function but evaluates directly the energies and transi- 



tion matrix elements for the transfer of one or more nucle- 
ons. Another important point is that the self-energy-the cen- 
tral component of the formalism-has been shown to be an ex- 
act optical potential. On the one hand, one can employ the 
formalism for pure ab initio studies. On the other, the strong 
link with the response to experimental probes can be used to 
constrain and improve phenomenological models. An exam- 
ple of this approach is the dispersive optical model recently 
derived for chain of calcium isotopes ll49il50ll . This global op- 
tical potential reproduces with high accuracy the known elas- 
tic scattering data, up to energies of 200 MeV. Thus, Green's 
functions hold a promise of both bridging nuclear structure 
and reactions and for connecting the (relatively few) isotopes 
amenable of ab initio calculations to the rest of the nuclear 
landscape. 

In practical applications of Green's functions theory, one 
expands the self-energy in terms of resummations of Feynman 
diagrams and truncates the series in a way that allows for fur- 
ther systematic improvements of the formalism. A powerful 
scheme for non-perturbative expansions is the so called Fad- 
deev random phase approximation (FRPA) that explicitly ac- 
counts for particle-vibration couplings ll28il5lll . First applica- 
tions of this approach were devoted to ie O. One single calcu- 
lation yielded the basic information to be used for microscopic 
studies of spectroscopic factors 1291152 ], the excitation spec- 
trum l30ll . two-nucleon emission l53l 54ll and the nucleon- 
nucleus optical potential f55ll . Thus, the FRPA method pur- 
sues a global description of the many-body dynamics, far be- 
yond ground state properties alone. Another difference with 
the coupled-cluster approach is that Green's functions allow, 
via a diagrammatic approach, to directly introduce the corre- 
lations outside the model space that are associated with short- 
range degrees of freedom I56ll57ll . It can therefore be applied 
to interactions with strong short-range cores (for example, the 
Argonne model f58ll was employed in Ref. ll3llo . 

Self-consistent FRPA calculations of 16 were first per- 
formed by the authors of Ref. 1I29I1 . and subsequently extended 
to fully ab initio calculations in spaces up to eight oscillator 
shells HH] . To our knowledge, the combination of the ran- 
dom phase approximation phonons and the proper treatment 
of the energy dependence of the interaction vertex, makes this 
the most accurate evaluation of single-particle states available 
for this nucleus. In this work we extend the range of applica- 
tions of this formalism to studies of quasiparticle states around 
56 Ni, as a first application of Faddeev random phase approx- 
imation to 1 pOf shell nuclei. As mentioned above, there is 
quite some experimental interest in single-particle properties 
around 56 Ni. 

This work is organized as follows. Sec. iHl re views the Fad- 
deev random phase approximation approach and discusses the 
approximations made to calculate the self-consistent propa- 
gator. The convergence of single-particle properties is dis- 
cussed in Sec. Hill and the results for the spectral function are 
described in Sec. [IV] We refer the reader directly interested to 
discussion of physics results to the latter section. Conclusions 
are drawn in Sec. [V] 
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H. FORMALISM 

This section serves as an overview of the formalism we em- 
ploy. The treatment of short-range physics and the implemen- 
tation of self-consistency, which are improved with respect to 
previous works, are discussed in details. The Faddeev random 
phase approximation expansion is also introduced, but we re- 
fer the reader for more details in Refs. [28, Sh- 
in the framework of Green's function theory, the object of 
interest is the single-particle propagator, instead of the many- 
body wave function. In the following, greek indices a,/3,. . . 
label the orthonormal basis set of single-particle states in- 
cluded in the model space, while latin indices k,l, , . . refer 
to many-body states. We employ the convention of summing 
over repeated indices, unless specified otherwise. The single- 
particle propagator can be written in the so-called Lehmann 
representation as l23l l59ll 
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where X" a = <^ +1 l4l^o> <X = <^ -1 |Ca,|^» are the spec- 
troscopic amplitudes, c a (c' a ) are the second quantization an- 
nihilation (creation) operators and e+ = E^ +l - (s7 = 
- Ef- 1 ). With these definitions, and IT^ 1 ) are the 

eigenstates, while E^ +l and E^~ l are the corresponding ener- 
gies of the (A + l)-nucleon system. Therefore, the poles of the 
single-particle propagator reflect the energy transfer observed 
in pickup and knockout reactions. 

The single-particle propagator g a p{w) enters the Dyson 
equation as 



(2) 



It depends on the irreducible self-energy 1,*{u). The latter can 
be written as the sum of two terms. The first terms describes 
the average mean-field (MF) while the second term contains 
dynamic correlations, 
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In Eq. g°(a>) is the so-called unperturbed single- 

propagator, corresponding to nucleons moving under the ef- 
fect of the kinetic energy part of the total Hamiltonian. The 
localization of the single-particle states in the nuclear mean- 
field is due to the term "E MF , which extends the Hartree-Fock 
potential to that of a fully correlated density matrix. The term 
V a p,ys represents the antisymmetrized matrix elements of the 
nucleon-nucleon interaction. In this work, these will be ap- 
proximated by an effective interaction, Eq. ©, discussed in 
Sec. Ill Al Eq. (0) introduces the two-particle-one-hole (2plli) 
and two-hole-one-particle (2hlp) irreducible propagator R(oj). 
In its Feynman expansion, this propagator contains all dia- 
grams with any number of particle and hole lines except those 
that allow the intermediate propagation of one single line. It 
can therefore be interpreted as carrying the complete infor- 
mation on all configurations that cannot be reduced to a nu- 
cleon interacting with the average nuclear field. In particular, 



it includes the coupling of single-particle states to collective 
vibrations like giant resonances. 

The theoretical spectroscopic factors and Z„ for removal 
and addition, respectively, of a nucleon, are given by the nor- 
malization integral of the corresponding overlap wave func- 
tions. In the notation of Eq. (Q]), these are 



a a 
a a 

The hole states are normalized according to 



(4) 



dZ* afS (a>) 



da> 



(5) 



which follows directly from the Dyson equation @. The same 
relation applies to particle states, with 3/* replaced by (<Y£)*. 
Because of the analytical properties of £*(«), the derivative 
term results always in a negative contribution, leading thereby 
to a quenching of the spectroscopic factors. 

It must be stressed that Eqs. (f2]) and ([3]) do not involve 
any approximation. Therefore, the full knowledge of 1, MF 
and R{a>) would be equivalent to the exact solution of the 
Schrodinger equation. In practical calculations, it is always 
necessary to truncate the full space of available single-particle 
states to a finite model space and to select a limited set of 
many-body correlations. The approximations employed in the 
present work to evaluate these quantities are discussed in the 
rest of this section. 



A. Short-range physics and effective interaction 

The present calculations were performed within a large 
but finite set of harmonic oscillator states, including single- 
particle states up to ten major shells. In order to treat the short- 
range part of the nucleon-nucleon interaction, one must resum 
explicitly the series of ladder diagrams for two nucleons out- 
side this model space. These contributions are included in 
the self-energy in two different ways. Firstly, they are explic- 
itly added to the mean-field part, £ MF , in order to reconstruct 
the contribution of short-range correlations in the full Hilbert 
space. Secondly, they are included in a regularized effective 
Hamiltonian which is used to calculate the long-range part 
of correlations — described by R(a>) — inside the chosen model 
space. This app roach leads to calculating the well known G- 
matrix J22i l60ll . which is then used as an energy dependent 
effective interaction inside the model space l8lll . In this case 
the Y. MF part of the self-energy will also depend on energy. 
This is given by 

yS 
yd k 
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where G a p ty §(u>) are the matrix elements of the G-matrix inter- 
action. Eq. (0 differs from the standard Brueckner-Hartree- 
Fock potential in the fact that the mean-field is not represented 
by a set of independent nucleons filling various orbits. Rather, 
the medium is described by the hole spectroscopic amplitudes 
J/* that lead to the fully correlated density matrix. Since the 
latter are obtained by solving the Dyson equation, they must 
be obtained iteratively in a self-consistent way. Whenever the 
second term in Eq. (0 is neglected, this procedure simply re- 
duces to solving the standard Brueckner-Hartree-Fock equa- 
tions. However, as soon as the response of the medium R(a>) 
is accounted for, the single-particle propagator becomes frag- 
mented and H MF describes the interaction of a particle with 
the "correlated" medium. 

We remind that using the G-matrix in Eq. (0 corresponds 
to summing the mean-field term form Eq. (0 and ladder di- 
agrams with intermediate two-particle states outside the cho- 
sen model space. This partitioning procedure has two main 
consequences. First, the effects of short-range physics at the 
two-body level on the total energy are included in the renor- 
malized interaction. This leads to a softer force that can be 
applied within a "low-momentum" model space. Secondly, 
due to the explicit energy dependence, the term Y, MF (u>) con- 
tributes as well to the normalization of spectroscopic factors, 
Eq. ([5]). This provides a natural way to determine the amount 
of strength that the free interaction would admix into config- 
urations outside the model space, see for example the discus- 
sion of Refs. I56l,l57ll . Thus, the present approach differs from 
methods based on the renormalization group, where instead 
an explicit renormalization of the effective operators would 
be required l6lll . 

The energy dependence of the G-matrix becomes cumber- 
some in calculating the polarization propagator. We define 
therefore a static effective interaction for our model space to 
be used in our calculations of the second term on the right- 
hand side of Eq. (0. To do this we evaluate the average energy 
for the harmonic oscillator single-particle states according to 
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where the sum is limited to those states that correspond 
to filled orbits in the independent particle model. Note 
that in Eq. (0 the single-particle energies are derived itera- 
tively while the oscillator wave functions remain unchanged. 
Clearly, these orbits are a crude approximation to the real 
quasiparticle states and will not be used to construct R(co) in 
Sec. Ill Bl -which will be rather determined in a self-consistent 
fashion. The main purpose of the above procedure is to yield a 
prescription for obtaining a starting energy independent effec- 
tive interaction. Following Gad and Miither 16211 . we use the 
single-particle energies from Eq. (0 to define an interaction 
for the given model space by 



Vpft* = -[G p, yS {cj = 4° + s'p + G ap , y5 {oj = 4° + 4°)] . 

__ (8) 

The G-matrix can be computed according Ref. 1221 16011 for 
negative energies, up to about -5 MeV. For larger values we 



fix the starting energy to a> — -5 MeV. Note that the starting 
energy appearing in Eq. (0 is shifted by the quasihole poles 
s7. This ensures that the energy dependence of 'L MF (co) is 
fully accounted for all quasiparticle states in the IpOf shell. 



B. The Faddeev random phase approximation method 

The polarization propagator R(to) can be expanded in terms 
of simpler Green's functions that involve the propagation of 
one or more quasiparticle states. This approach has the ad- 
vantage that it aids in identifying key physics ingredients of 
the many-body dynamics. By truncating the expansion to a 
particular subsets of diagrams or many-body correlations, one 
can then construct suitable approximations to the self-energy. 
Moreover, since infinite sets of linked diagrams are summed 
the approach is non-perturbative and satisfies the extensivity 
condition fTill . This expansion also serves as a guideline for 
systematic improvements of the method. 

Following Refs. j28ll5lll . we first consider the particle-hole 
polarization propagator that describes excited states in the A- 
particle system 
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and the two-particle propagator that describes the addi- 
tion/removal of two particles 



5 afi,yS 



(CO) = J] 



jN+2 



c y c S^0 ) 



(E™-E») + k, 



y (X\cW s \^- 2 ) (^- 2 \ C/3 c a \X) 

Z—i (rpN I?N~2\ 

k *> - [E Q -E k j-nj 



(10) 



These Green's functions contain in their Lehmann represen- 
tations all the relevant information regarding the excitation of 
particle-hole and two-particle or two-hole collective modes. 
In this work we are interested in studying the influence of gi- 
ant resonance vibrations, which can be described within the 
random phase approximation (RPA). In the Faddeev RPA ap- 
proach, the propagators of Eqs. (0 and ( fTOb are then evalu- 
ated by solving the usual RPA equations, which are depicted 
diagrammatically in Fig.Q] Since these equations reflect two- 
body correlations, they still have to be coupled to an addi- 
tional single-particle propagator, as in Fig.0 to obtain the cor- 
responding approximation for the two-particle-one-hole and 
two-hole-one-particle components of R(oj). This is achieved 
by solving two separate sets of Faddeev equations, as dis- 
cussed in Ref. j28ll . 

Taking the two-particle-one-hole (2plh) case as an exam- 
ple, one can split R ( - 2pl,i) (to) in three different components 
R('\u)) (i = 1, 2, 3) that differ from each other by the last pair 
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FIG. 1 : (Color online) Diagrammatic equations for the polarization 
(above) and the two-particle (below) propagators in the RPA ap- 
proach. Dashed lines are matrix elements of the effective nucleon- 
nucleon interaction, Eq. l[8j. The full lines represent the independent- 
particle model propagator g IPM (aj), which is employed instead of the 
fully dressed one. See the text for details. 



of lines that interact in their diagrammatic expansion, 

Ra/ly^v^ = [ G °«/iy,//vi(w) - G°^ ayjtyX ((J)]+ R%y, M vA^ ' 

1=1,2,3 

(11) 

where G° (u>) is the 2plh propagator for three freely propa- 
gating lines. These components are solutions of the following 
set of Faddeev equations Ir33ll 

x + tfbwM ( 12 ) 

+ G ^, v ,, A ,,^ vA (a))-G v ,, f j,, A ,, llvA (u>)\, i =1,2,3 

where (i,j,k) are cyclic permutations of (1,2,3). The in- 
teraction vertices r w (w) contain the couplings of a particle- 
hole (ph), see Eq. (0, or two-particle/two-hole (pp/hh), see 
Eq. ([Tol l, collective excitations and a freely propagating line. 
The propagator R(co) which we employ in Eq. (|3]l is finally 
obtained by 

where the matrix U has the effect of renormalizing the 
strength of the dynamic self-energy. This correction ensures 
consistency with perturbation theory up to third order. The 
explicit formulae of the matrices T w (w) and U are given in 
terms of the propagators of Eqs. (0, ( fTOb and (14\ . and the 
interaction V a ^ y s. They are discussed in detail in Ref. l5lll . 
The calculation of the 2hlp component of R(co) follows com- 
pletely analogous steps. 

The present formalism includes the effects of ph and pp/hh 
motion simultaneously, while allowing interferences between 
these modes. These excitations are evaluated here at the 
RPA level and are then coupled to each other by solving 
Eqs. ([T2l i. This generates diagrams as the one displayed in 
Fig. |2] The Faddeev equations also ensure that the Pauli prin- 
ciple is correctly taken into account at the 2plh and 2hlp 




FIG. 2: (Color online) Example of one of the diagrams that are 
summed to all orders by means of the Faddeev random phase approx- 
imation Eqs. d 12b (left). The corresponding contribution to the self- 
energy, obtained upon insertion into Eq. (O, is also shown (right). 



level. In addition, one can in principle employ dressed single- 
particle propagators in these equations to generate a fully self- 
consistent solution, as done in Refs. l29[|3lll for valence orbits 
around 16 0. 



C. Self-consistent approach 

In the self-consistent Green's function approach, both the 
I. MF part of the self-energy and the polarization propaga- 
tor R(cl>) are expressed directly in terms of the exact single- 
particle propagator g(a>). The lines in Figs. Q] and [2] should 
thus represent the fully dressed propagator obtained by solv- 
ing the Dyson equation. Since the degrees of freedom con- 
tained in Eq. (Q]i are excitations of the fully correlated sys- 
tem, the formalism does not depend on an explicit reference 
state. Normally, one first computes Eq. <j3j in terms of an ap- 
proximate propagator. The solution of Eq. (O is then used to 
calculate an improved self-energy and the procedure is iter- 
ated to convergence. Baym and Kadanoffhave shown that the 
self-consistency requirement implies the conservation of both 
microscopic and macroscopic properties 164 16511 . Intuitively, 
the self-consistency requirement becomes important when- 
ever dynamical correlations modify substantially the response 
with respect to the Hartree-Fock mean-field (an example is 
the band-gap error problem in diamond crystals l66ll ). When 
applying standard Hartree-Fock theory to nuclear structure, 
most realistic interactions predict unbound nuclei and valence 
orbits in the continuum. This is a very poor starting point 
for any application of perturbation theory and other many- 
body techniques. However, the self-consistent approach re- 
quires using correlated quasiparticle energies and wave func- 
tions [the poles and residues of Eq. ([TJ]. These degrees of 
freedom form an optimal starting point for studies of many- 
body dynamics at the Fermi surface. 

Accounting for the fragmentation of the single-particle 
propagator in the Faddeev random phase approximation in- 
creases the computational load as one moves to larger nuclei 
and model spaces. In this situation it is convenient to expand 
R(u)) in terms of an independent-particle model (IPM) prop- 
agator. This should approximate the dressed one but with a 
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limited number of poles. Thus, we solve Eq. ( fT2t in terms of 
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where F represent the set of occupied orbits. The single- 
particle energies s 1PM and wave functions <f> are chosen such 
that g IPM (u>) coincides with the real propagator g(a>) at the 
Fermi surface. To do this we define the following moments of 
the poles of Eq. ([TJ. 



M „ y (gr g v jg (vff 



(15) 

where £> is the Fermi energy. Eq. (fT4b is determined by im- 
posing M^ PM = and M^' PM = The purpose of 
Eq. <fT3T > is to define a set of effective single-particle orbits 
and energies that conserve the total spectroscopic strength car- 
ried by the self-consistent propagator and the centroids of its 
fragmented states. While effective single-particle properties 
form an appropriate starting point to evaluate R(oj), it remains 
clear that they only represent average quantities. Instead, it is 
Eq. ([T]i that must be related to experiment. 

The propagator g ,PM (a>) is derived from Eq. ([TJ and it still 
needs to be evaluated in a iterative way. Therefore, the result- 
ing propagator R(cj) is (partially) self-consistent. We stress 
that Eq. © can be calculated easily from the fully dressed 
propagator. Thus self-consistency is achieved exactly at the 
mean-field level. 



III. CALCULATIONS AND CONVERGENCE 

The present calculations were performed using a harmonic 
oscillator basis and including up to ten major harmonic os- 
cillator shells. We label these spaces with N max - 3, 5, 7, 
or 9, where N - 2n + I. For the largest model space em- 
ployed, Nmax - 9, all partial waves with orbital angular mo- 
mentum I < 7 were included. This amount to 368 single- 
particle states for each particle species, protons and neutrons 
in our case. The total number of available Slater determinants 
for 56 Ni without any particular restrictions is proportional to 
the product of the two binomials 



368 
28 



368 
28 



a number which clearly exceeds the capabilities of any direct 
large-scale diagonalization procedure. 

The codes utilize a y'y'-coupling scheme to decouple the 
Faddeev equations of Eq. (fT2l . At each iteration, the RPA 
equations are solved in the particle-hole and the particle- 
particle and hole-hole channels using the single-particle or- 
bits and energies from Eq. ( fT~4T >. The resulting propagators 
are inserted in Eqs. (ITZb . which in turn are casted into a 
non Hermitian eigenvalue problem J28|. In our largest cal- 
culation we diagonalize dense matrices of dimension up to 



26 x 10 3 two-particle-one-hole states. This number is bound 
to increase when larger nuclei are investigated or more details 
of nuclear fragmentation (that is more poles) are included in 
g IPM (u>) (fT4l i. The numerical implementation of the Faddeev 
random phase approximation required careful optimization in 
evaluating the elements of the Faddeev matrix and a proper 
generalization of the Arnoldi algorithm l67ll to employ mul- 
tiple pivots. Similar improvements allowed to extend large 
scale calculations from the A — 16 mass region to the A = 56 
mass region. The dimensions reached in this work represent 
roughly the upper limit when using table top single processor 
computers. Obviously, there is much to gain by taking advan- 
tage of modern supercomputer facilities and future research 
efforts should be put into parallelization of the present algo- 
rithms. 

Model spaces of eight to ten major shells are large enough 
for a proper description of the response due to long-range cor- 
relations. These include excitations of several MeVs into the 
region of giant resonances. The effects of short-range physics 
are also included by using a G-matrix and an effective interac- 
tion as discussed in Sec. Ill Al These are derived using the chi- 
ral nucleon-nucleon interaction N 3 LO by Entem and Mach- 
leidt ll47ll . This interaction employs a cutoff of A = 500 MeV. 

Typical realistic two-nucleon interactions fail in reproduc- 
ing the spin-orbit splittings and gaps between different shells. 
In particular, for the N - 2% and Z = 28 subshell closures 
these lead to an underestimation of the gap at the Fermi sur- 
face. In these cases, a complete diagonalization of the Hamil- 
tonian would predict a deformed ground state of these nuclei 
even when they are experimentally known to be good spher- 
ical closed-shells systems ifTTl |68[| . This issue can be cured 
with a simple modification of the monopole strengths of the 
interaction. Recently, Zuker has reported that the same cor- 
rection works well for several isotopes throughout the nuclear 
chart and proposed that this may be interpreted as a signature 
of missing three-nucleon interactions II69I1 . 

The inclusion of three-nucleon interactions to the Faddeev 
random phase approximation formalism is beyond the scope 
of the present work. However, it will be shown in Sec. I III Al 
that properly reproducing the Fermi gap is crucial in order to 
obtain meaningful results for the valence space spectroscopic 
factors. Thus, we follow Ref. l69ll and modify the monopoles 
in the N 3 LO interaction model as 



AVj r 
AV T f 



AV T fr -(-l) T K M , 
AVf f — 1.5(1 — T)km 



' ff ~ T L ^ y Si ~ y ~'^ L ~ ± >* M > (16) 
where / and r stand for the Brueckner-Hartree-Fock states as- 
sociated to the O/7/2 and the (Ip3/2,lpi/2,0f5/z) orbits, respec- 
tively. In the limit of large spaces, the Brueckner-Hartree- 
Fock orbits converge to the Hartree-Fock states and this cor- 
rection becomes independent of the choice of the single- 
particle basis. Note that the prescription of Eq. ( fT&b modifies 
only a few crucial matrix elements while about six millions of 
them are defined in the N max = 9 model space. 

We also note that the present truncation of the model space, 
in terms of the number of oscillator shells, does not separate 
exactly the center of mass motion. Coupled-cluster calcula- 
tions have shown that the error introduced by this truncation 
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FIG. 3: (Color online) Dependence of neutron spectroscopic factors 
(given as a fraction of the independent-particle model value) for the 
I Pi 12 and the O/7/2 valence orbits with respect the ph gap AE p i,. For 
each model space, different points correspond to different choices of 
Km in the range 0.4 - 0.7 MeV. 



becomes negligibly small for large model spaces such as the 
ones employed here and therefore it does not represent a major 
issue jr7ll70ll . In calculations of binding energies, it custom- 
ary to subtract the operator for the kinetic energy of the center 
of mass directly from the Hamiltonian. This term automati- 
cally corrects for the zero point motion in oscillator basis but 
it depends explicitly on the number of particles. In this work, 
we are interested in transitions to states with different num- 
bers of nucleons [(A+l] and aim at computing directly the 
differences between the total energies. Therefore, the above 
correction should not be employed in the present case. One 
may note that the separation of the center-of-mass motion is 
an issue related to the choice made for the model space, rather 
than the many-body method itself. For example, expressing 
the propagators directly in momentum space would allow an 
exact separation. In this situation, the transformation between 
the center-of-mass and laboratory frames for systems with a 
nucleon plus a A nucleons (or (A - I) nucleons) core would 
also be simple. 



A. Choice of k m 

Eq. ( [Tol l introduces a single parameter (km) in our calcu- 
lations. The reason for this modification is that the spectro- 
scopic factors of the valence orbits are strongly sensitive to 
the particle-hole gap. This sensitivity is to be expected since 
collective modes in the 56 Ni core are dominated by excitations 
across the Fermi surface. Smaller gaps imply lower excitation 
energies and higher probability of admixture with valence or- 
bits. In order to extract meaningful predictions for spectro- 
scopic factors it is therefore necessary to constrain the Fermi 
gaps for protons and neutrons to their experimental values. 

To investigate this dependency we repeated our calcula- 
tions for values of km in the range 0.4 - 0.7 MeV. Fig. [3] 
shows the resulting neutron spectroscopic factors for the va- 
lence P3/2 quasiparticle and fi/2 quasihole. These are plot- 
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FIG. 4: (Color online) Dependence of the neutron l/?3/2 particle en- 
ergy and the O/7/2 hole energy with respect to the oscillator frequency 
and the size of the model space. 



ted as a function of the calculated particle-hole gap AE p h = 
e U>3n n=o ~ s Qfy2 k=o- The results correspond to model spaces 
of different dimensions (eight or ten oscillator shells) and os- 
cillator frequencies (hQ. = 10 or 18 MeV). The gap AE p h in- 
creases with km but the dependence on the model space is 
weak. We notice that, once the experimental value of AE p i, 
is reproduced, the spectroscopic factors are well defined and 
found to be converged with respect to the given model space. 

All results reported below were obtained with a fixed value 
of k m = 0.57 MeV. In the N max = 9 model space and an 
oscillator energy hQ. = 10 MeV, this choice reproduces the 
experimental gaps at the Fermi surface for both protons and 
neutrons to an error within 70 keV. From Fig. [3] one infers 
that the calculated spectroscopic factors are reliable to within 
1 - 2% of the independent-particle model value. 



B. 



Convergence with respect to choice of model space 



Fig. [4] shows the dependence of the neutron lpj,/2 particle 
and the O/7/2 hole energies with respect to the oscillator fre- 
quency and the size of the model space. As can be seen from 
this figure, the single-particle energies for these two single- 
particle states tend to stabilize around eight to ten major shells. 
This finding concords both with coupled-cluster calculations 
that employ a G-matrix as effective interaction for 16 0, see 
Refs. |70l UM > an d with analogous Green's functions stud- 
ies 13111 - It remains however to make an extensive comparison 
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hQ [MeV] 

FIG. 5: (Color online) Dependence of neutron single-particle ener- 
gies on the oscillator frequency. The energies plotted here are the 
poles of g(oS) corresponding to the valence 1 pOf orbits. Calculations 
were performed for N„ wt = 9 and I < 7. 



between coupled-cluster theory and the Green's function ap- 
proach in order to find an optimal size of the model space with 
a given nucleon-nucleon interaction. Finally, we plot in Fig. [5] 
the neutron valence single-particle energies for all the single- 
particle states in l/?0/-shell. The latter results were obtained 
with our largest model space, ten major shells with N max — 9 
and the single-particle orbital momentum I < 7. As can be 
seen from this figure, there is still, although weak, a depen- 
dence upon the oscillator parameter. To perform calculations 
beyond ten major shells will require non trivial extensions of 
our codes. 



IV. RESULTS FOR THE SPECTRAL FUNCTION 

Our results for spectroscopic factors for the 1 pO/-shell va- 
lence orbits and the corresponding single-particle energies are 
collected in Table Q] In general, the modified N 3 LO interac- 
tion predicts single-particle energies about 2-3 MeV lower 
than the experimental ones. The Coulomb shift between cor- 
responding neutron and proton orbits is calculated to be about 
10.2 MeV and it is closer to the empirical value of 9.5 MeV. 
For the oscillator parameter chosen, HQ. = 10 MeV, we obtain 
an inversion of the lpi/2 and O/5/2 excited states in 57 Ni, with 
respect to the experiment. However, this discrepancy disap- 
pears for larger values of ftQ (see Fig. [5J. This effect is in 
fact smaller than the residual dependence on the model space 
and therefore no conclusion can be made about the ordering 
for the fully converged result. The spectroscopic factor for the 
transition between the ground states of 57 Ni and 56 Ni was ex- 
tracted from high energy knockout reactions in Ref. J2]. The 
self-consistent Faddeev random phase approximation result 
for this quantity yield 65% of the independent-particle model 
value, and agrees with the empirical data within experimen- 
tal uncertainties. The theoretical spectroscopic factors for the 
excited states in 57 Ni are similar, with the O/5/2 state being 



s + n ,e' k Z„/(2j+\),Z k /(2j+l) 





FRPA 


Exp. 


FRPA 


57 Ni: 








vl.Pl/2 


-11.43 


-9.134 


0.63 


vO/s/2 


-10.80 


-9.478 


0.59 


vlp 3/2 


-12.78 


-10.247 


0.65 


55 Ni: 








v0/ 7/2 


-19.22 


-16.641 


0.72 


57 Cu: 








nlpi/2 


-1.28 


+0.417 


0.66 




-0.58 




0.60 


nlp 3/2 


-2.54 


-0.695 


0.67 


55 Co: 








TO/7/2 


-9.08 


-7.165 


0.73 



Exp. 



TABLE I: Energies (in MeV) and spectroscopic factors (as a fraction 
of the independent-particle model) for transitions to the 1 pOf valence 
orbits, obtained for TiQ = 10 MeV, N„ mx = 9 and k m = 0.57 MeV. 
The experimental single-particle energies are taken from |72Tl . The 
measured spectroscopic factor for transfer between the ground states 
of 57 Ni and 56 Ni is from Ref. 0]. 



somewhat smaller, at about 59%. A larger value is obtained 
for knockout to the ground state of 55 Ni, which is predicted 
to be 72%. The results for proton transfer to particle (hole) 
states in 57 Cu ( 55 Co) are only slightly larger. According to the 
analysis of Fig. [3] it is expected that these predictions are con- 
verged within 1-2% of the independent-particle model values. 

Past studies H HI 111] have questioned whether low-energy 
quasiparticle states in 57 Ni are strongly admixed to excitations 
of a soft 56 Ni core. The results obtained here do not sug- 
gest substantial differences with respect to other known closed 
shell nuclei. The spectroscopic factors from Tableware in line 
with observations from stable nuclei J^l p74Tl and support the 
hypothesis that 56 Ni is a good closed-shell nucleus. In our 
calculations we find that the l/?3/2 quasiparticle state of 57 Ni 
carries 65% of the strength for this orbit. Another 20% is 
located in the particle region below MeV (above this 

energy strength associated with the 2p\f0h shell starts to ap- 
pear), and about 3% is in the hole region above s k = -40 
MeV (see Fig. [6j. Similarly, the O/7/2 state has 72% of the 
independent-particle model strength in the quasihole peak (the 
ground state of 55 Ni), 10% in the fragmented hole region, and 
3% in the fragmented particle region. This analysis confirms 
that the main mechanism responsible for the quenching of the 
spectroscopic factors lies in the admixture between single- 
particle states and collective excitations in the region of giant 
resonances [24]. Due to these correlations a large part of the 
missing strength from the valence peak is shifted and spread 
over an adjacent region about 15 - 20 MeV wide. Further 
reduction of the spectroscopic factors comes from the mixing 
with configurations at much higher energies and momenta and 
is accounted for through the energy dependence of Eq. (|5J, 
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FIG. 6: (Color online) Spectral strengths for one-neutron transfer on 56 Ni obtained from the self-consistent single-particle propagator g(w). 
Poles above (below) the Fermi energy, E F , correspond to transition to eigenstates of 57 Ni ( 55 Ni). The respective spectroscopic factors are 
given as a fraction of the independent-particle model value. The quasiparticle poles corresponding to the valence orbits of the 1 pOf shell are 
indicated by arrows. A logarithmic scale was chosen to put in stronger evidence the distribution of the fragmented strength. Results are for 
nn=W MeV, N max =9 and k„=Q.51 MeV. 



The information carried by the calculated single-particle 
propagator g(u>) is collected in Fig.|6]for neutrons and partial 
waves up to 1=3. The plot shows the spectral strength asso- 
ciated with each pole of g(aj). Fragments below the Fermi 
surface (£» refer to the separation of a neutron ( 55 Ni), while 
those above correspond to neutron addition ( 57 Ni). The poles 
corresponding to the IpOf valence orbits are indicated by ar- 
rows. These are the same single-particle energies that have 
been discussed in Sec. 1111 B I and Tab. U As noted above, 
the fragments found at slightly higher energies (just above 
~ -8 MeV) originate from the mixing of these orbits 
with two-particle-one-hole configurations and collective ex- 
citations of the nucleus. The overall fragmentation effect is 
substantial but not strong enough to destroy the single-particle 
character of the principal quasiparticle peaks (note that a log- 
arithmic scale has been chosen in Fig. [6] in order to make 
smaller fragments of the spectral distribution more visible). 
This observation supports the use of valence single-particle 
states as the relevant degrees of freedom that govern low-lying 
excitations, as assumed in conventional shell-model applica- 
tions. We note that the question whether a system can be 
approximated as a good shell closure is better addressed by 



analyzing the spectroscopic factors and strength distribution 
rather than occupation numbers, since the latter are integrated 
quantities l82il . While unoccupied states can be probed by 
the addition of a nucleon, occupied states are accessed by 
knockout to states of the A - 1 nucleon system. A similar 
fragmentation pattern is therefore seen for the O/7/2 orbit but 
reversed below the Fermi surface. Interestingly enough, the 
Faddeev random phase approximation predicts that states cor- 
responding to orbits in the IsOd and 2 si dOg shells maintain 
a strong single-particle character even though they are further 
apart from the Fermi surface. The fragmentation of these or- 
bits requires excitations across shells of different parity (e.g. 
1 sOd and 1 pOf) and could become stronger if the energy dif- 
ference among major shells is reduced. Indeed, a compari- 
son of our results with electron scattering measurements on 
58 Ni lF75h suggests that the N 3 LO interaction tends to over- 
estimate the gaps between these major shells. Note that in 
the present calculations the 2sld0g quasiparticles are found 
at energies of about -3 MeV and overlap with the fragmented 
IpOf states. 

Far from the Fermi energy Ep, the mixing with complex 
configurations becomes strong and it is no longer possible to 
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identify sharp quasiparticle and quasihole states. Still, the 
energy region occupied by the major shells can be identi- 
fied clearly. The N 3 LO interaction places the states associ- 
ated with the IsOd shell between -60 and -30 MeV, while 
the 0/7-shell states appear below -50 MeV. Other hole frag- 
ments are observed around -30 MeV for the lpi/2, l/?3/2 and 
O/5/2 partial waves. These originate from particle states in 
the 1 pOf shell that are partially occupied due to the smear- 
ing of the Fermi surface. Nucleon knockout from these orbits 
requires little energy transfer and leads to low-lying states in 
55 Ni or 55 Cu. These states originate from the mixing of par- 
ticle orbits with two-hole-one-particle configurations (Bp tA P>') 
in the Dyson equation. Still, the 55 Ni ground state is strongly 
influenced by the O/7/2 hole component. 

Analogous fragmentation patterns extend to the shells fur- 
ther away from the Fermi surface, although these are not 
shown in Fig. [6] On the particle side, e+ = MeV marks 
the threshold for the single -particle continuum in the (A + 1)- 
nucleon system. Above this, the exact spectral function be- 
comes a continuous function of energy. In the present calcu- 
lations a structure of separate poles is found due to the dis- 
cretization of the model space. A continuum distribution also 
develops for the hole part of the exact spectral function below 
the energy s^=(Eq - Efzl,) - where S^~ l is the one- 

nucleon separation energy from the ground state of A - 1 par- 
ticles. The distribution on both sides of the Fermi surface is 
similar but not fully symmetric, the strength being stronger at 
large positive energies. This is because the (A + l)-nucleon 
system can access a larger phase space than a single hole 
within the A-nucleon ground state. This asymmetry is already 
observed at the level of the self-energy T.*(u>), a result in line 
with available fits of global optical potentials ll49ll50ll . 

For the proton case, the poles of g(a>) correspond to the ad- 
dition (removal) of a proton to the eigenstates of 57 Co ( 55 Cu). 
The corresponding spectral strength is substantially the same 
as that discussed for neutrons due to the almost exact isospin 
symmetry of the nuclear force. However, it is shifted to higher 
energies by the Coulomb repulsion. 



V. CONCLUSIONS 

The aim of this work has been to extend large-scale calcu- 
lations of self-consistent Green's functions to medium mass 
nuclei and to investigate the properties of the single-particle 
spectral function of 56 Ni. 

Many-body Green's functions hold a number of interest- 
ing mathematical properties. Since one aims at obtaining 
excitations relative to a reference nucleus calculations scale 
more gently when increasing the number of particles as op- 
posed to direct large-scale diagonalization methods. Only 
connected diagrams are summed to all orders so that the 
extensivity condition is satisfied lfl4ll . Moreover, the self- 
consistent approach provides a path to ensure the conserva- 
tion of basic macroscopic quantities. However, the greatest 
advantage of the self-consistent Green's function formalism 
is that its building blocks, the many-body propagators, con- 
tain information on the response to several particle transfer 



and excitation processes. Therefore, they can be directly com- 
pared to a large body of experimental data. Due to these char- 
acteristics the formalism can be used to gain unmatched in- 
sights into the many-body dynamics of quantum mechanical 
systems. Within this framework, the Faddeev random phase 
approximation method proposed in Ref. |28| is a good candi- 
date to pursue ab initio studies of medium mass isotopes. 

In this work we have presented the basic details of calcu- 
lating the Faddeev random phase approximation expansion 
and discussed results for the spectral function of 56 Ni. This 
is the first application of the self-consistent Green's func- 
tion approach to the 1 pOf shell region. The calculations em- 
ploy the chiral N 3 LO two-nucleon interaction, with a modified 
monopole to account for missing many-nucleon forces. In ad- 
dition to this one-parameter modification of the Hamiltonian, 
the only remaining parameters that enter our calculations are 
those defining the nucleon-nucleon interaction. 

Calculations were performed in models spaces including up 
to ten major oscillator shells. These large spaces are large 
enough to allow for a sophisticated treatment of long-range 
correlations. The quasiparticle and quasihole energies of the 
1 pOf valence orbits were found to be rather well converged. 
In the largest calculations they appeared to be almost con- 
stant for oscillator frequencies in the range fiQ. e [8, 20] MeV. 
These convergence properties are possible thanks to a predi- 
agonalization of the effects of short-range correlations. This 
is done using the G-matrix technique ll22ll60Tl to resum ladder 
diagrams outside the model space. Our results put in evidence 
the strong sensitivity of spectroscopic factors on the particle- 
hole gap at the Fermi surface. For the N = 28 and Z = 28 
subshell closures the bare N 3 LO potential fails in describing 
the experimental gap (in an analogous way to other realistic 
two-nucleon interactions ifTTIl ). This effect has been attributed 
to missing three-nucleon interactions J69ll . It is found that a 
proper correction of few monopole terms of the Hamiltonian 
allows us to extract reliable results for the fragmentation of 
single-particle strength. 

Fully self-consistent Faddeev random phase approximation 
calculations have till now only been presented for 16 0. The 
extension to accurate ab initio calculations in the 1 pOf shell 
represent a major technical advance. However, no substantial 
use of parallel computation has been made in applying this 
formalism. Improvements in numerical algorithms are still 
possible and it is expected that they will allow a better treat- 
ment of fragmentation in the self-consistent approach, as well 
as pushing the limits of present calculations well beyond mass 
A = 56. Another obvious extension is the inclusion of explicit 
three-nucleon forces. Within the framework of self-consistent 
Green's function theory this has already been achieved for nu- 
clear matter studies 117611 . Similar developments can be ex- 
pected for finite nuclei as well. 

For open shell systems with weakly bound states and/or 
resonances, one needs a single-particle basis which can han- 
dle continuum states, as done in ll77ll . Normally, this leads 
to a much larger space and may require parallelization of 
our codes. On the other hand, the effective interaction 
among valence-space quasiparticles is already generated in 
the present calculations and can be used for standard shell- 
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model calculations in one or two major shells. The issue of 
degenerate unperturbed states for open shells systems has also 
been addressed for Green's functions theory in Ref. lf78ll by 
using a Bogoliubov-type quasi-particle transformation. In a 
self-consistent treatment, one may improve on this approach 
by extending the Faddeev-RPA method to include explicit 
configuration mixing between the nucleons inside the open 
shell. 

The N - 2d, and Z = 28 subshell closure has also at- 
tracted recent experimental interest following the discussion 
of whether the low-lying states of 57 Cu are strongly frag- 
mented due to a soft 56 Ni core, see for example Ref. [6fl . While 
no direct experimental information is available for the transi- 
tion between these two isotopes, the spectroscopic factor for 
neutron knockout from 57 Ni has been measured in Ref. 10]. 
The present calculations describe well the quenching of the 
experimental cross section. At the same time, we report pre- 
dictions for both proton and neutron transfer to the other va- 
lence orbits around 56 Ni. These calculations can thereby pro- 
vide theoretical benchmarks for the forthcoming experiments 



of Ref s. lHHH]. These spectroscopic factors are all in the range 
of 60%-70% of the independent-particle model value and in 
fair agreement with the observation of valence states in sev- 
eral stable nuclei j9]. The fragmentation pattern of valence 
orbits predicted by the Faddeev random phase approximation 
is also found in substantial agreement with what is known for 
closed shell nuclei l24ll and supports the description of 56 Ni 
as a doubly magic nucleus. Finally, we note that the effects of 
admixing configurations with several particle-hole excitations 
are not included in the present study. These effects can be ac- 
counted for by using configuration interaction (shell-model) 
methods. However, based on the analysis of Refs. 1 6 8, 179118011 
these corrections are not expected to be dominant. 
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